Quasi-two-dimensional optomechanical crystals with a complete phononic bandgap 
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A fully planar two-dimensional optomechanical crystal formed in a silicon microchip is used to create a 
structure devoid of phonons in the GHz frequency range. A nanoscale photonic crystal cavity is placed inside 
the phononic bandgap crystal in order to probe the properties of the localized acoustic modes. By studying the 
trends in mechanical damping, mode density, and optomechanical coupling strength of the acoustic resonances 
over an array of structures with varying geometric properties, clear evidence of a complete phononic bandgap is 
shown. 



Interest in cavity-optomechanical systems, in which light 
is used to sensitively measure and manipulate the motion of 
a mechanically compliant optical cavity, has grown rapidly 
in the last few years due to the demonstration of micro- and 
nano-scale systems in which the radiation pressure force of 
light is manifest. The physics of these systems is similar in 
many regards to the inelastic scattering of light by localized 
molecular vibrations (Raman scattering), where the mechan- 
ical resonance is now associated with a moveable end-mirror 
affixed to a spring (or hung as a pendulum), as is the case in 
a Fabry-Perot cavity. More complex cavity-optomechanical 
geometries, such as whispering-gallery mode structures | l]- 
[6l, nanomembranes placed within Fabry-Perot cavities, and 
near-field optical and microwave devices utilizing the gradient 
force 13 0, have also been demonstrated to produce strong 
radiation pressure effects. Two structures with particular en- 
gineerability of optical and mechanical properties are the pho- 
tonic crystal fiber (PCF) mHOl and the recently demonstrated 
chip-scale optomechanical crystal (OMC) 1 11 1. In the case of 
the OMC cavity, a combined photonic and phononic crystal 
formed in a thin nanobeam of silicon (Si) is used to create 
a 200 THz optical cavity simultaneous with a GHz acoustic 
cavity. 

Photonic crystals are engineered periodic dielectric struc- 
tures in which optical waves with wavelengths close to the 
period of the structure encounter strong dispersion, and in 
some cases are completely forbidden from propagating within 
a photonic bandgap frequency window. In the last decade. 
Si photonics has rapidly developed, in no small part due to 
the advent of silicon-on-insulator wafer technology in which 
a thin Si device layer sits atop a low-index insulating ox- 
ide layer. Planar Si photonic crystal circuits may be formed 
in the top Si device layer, which along with the integrabil- 
ity with micro-electronics, has provided an attractive setting 
for controlling photons. In a similar way, the periodic pat- 
terning of elastic structures can be used to control acoustic 
wave propagation and create phononic bandgaps. These struc- 
tures, known as phononic crystals, have been created in a va- 
riety of materials and over a wide range of scales and geome- 
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FIG. 1: a, Real space crystal lattice of the cross crystal with lat- 
tice constant a, cross length h, cross width w, and membrane thick- 
ness d. The bridge width is defined SiS b = a — h. b, Reciprocal 
lattice of the first Brillouin zone for the cross crystal, c, Phononic 
band diagram for the nominal cross structure with a = 1.265 ^um, 
h = l .220 jUm, w = 340 nm. Dark blue lines represent the bands with 
even vector symmetry for reflections about ihQx — y plane, while the 
red lines are the flexural modes which have odd vector mirror sym- 
metry about the x — y plane, d, Tuning of the bandgap with bridge 
width, b. Light grey, dark grey, and white areas indicate regions 
of a symmetry-dependent (i.e., for modes of only one symmetry) 
bandgap, no bandgap, and full bandgap for all acoustic modes, re- 
spectively. 
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tries IIT2HT91 for applications as diverse as the filtering and 
focusing of sound| 20 |, the processing of RF/microwave sig- 
nals |[2TH24I. and the earthquake proofing of buildings (251. 
Thin-film OMCs bring together photonic and phononic cir- 
cuitry 1 26- - 30 L enabling a new chip-scale platform for de- 
laying, storing, and processing optical and acoustic excita- 
tions f3T, "321. Building on our previous work in quasi- 
ID OMC systems |11|, we demonstrate here a quasi-two- 
dimensional OMC architecture capable of routing photons 
and phonons around the full two-dimensional (2D) plane of 
a Si chip, and enabling complete localization of phonons via 
a three-dimensional acoustic- wave bandgap. 

The phononic crystal used in this work is the recently pro- 



2 




+++++++++++++++++++ 
+++++++++++++++++++ 



+++++::::::::::::::::::::::::::::>+++ 
t+++H:':':*:':'x^^^^^^ 

++++++++++++++ 
+++++++++++4 ++ 
++++++++++++++ 
++++++++++++++^ 





1|Lim, 



dooooooooooooooooooooooooooo 
oooooooooooooooooooooooooooo 
ooooooooooooooooooooooooooo 
oooooooooooooooooooooooooooooo 
ooooooooooooooooooooooooooooo 
oooooooooooo oo o«o 0*0 000000000000 
• • • • • 

ooooooooooooo«o OlO 00 oooooooooooo 
ooooooooooooooooooooooooooooo 

oooooooooooooooooooooooooooooo 
ooooooooooooooooooooooooooooo 

oooooooooooooooooooooooooooooo 
ooooooooooooooooooooooooooooo 










1.5x10^ 



,00 s 



99 



97 



-20 -10 10 20 

3X (pm) 



o 
=3 

c 



FIG. 2: a, Scanning electron micrograph (SEM) of one of the fabricated 2D-0MC structures. The photonic nanocavity region is shown in false 
green color. In b, Zoom-in SEM image of the cross crystal phononic bandgap structure, c, Zoom-in SEM image of the optical nanocavity within 
embedded in the phononic bandgap crystal. Darker (lighter) false colors represents larger (smaller) lattice constant in the optical cavity defect 
region, d, FEM simulation of Ey electrical field for the optical cavity, e, Typical measured transmission spectra for the optical nanocavity, 
showing a bare optical g-factor of 2/ = 1.5 x 10^. 



posed 1291 "cross" structure shown in Fig. [T^. Geometri- 
cally, the structure consists of an array of squares connected 
to each other by thin bridges, or equivalently, a square lattice 
of cross-shaped holes. The phononic bandgap in this struc- 
ture arises from the frequency separation between higher fre- 
quency tight-binding bands, which have similar frequencies to 
the resonances of the individual squares, and lower frequency 
effective-medium bands with frequencies strongly dependent 
on the width of the connecting bridges, b = a — h \29\. A 
typical band diagram for a nominal structure (a = 1.265 /im, 
h = 1.220 /im, w = 340 nm) is shown in Fig.[T]:. Blue (red) 
lines represent bands with even (odd) vector symmetry for re- 
flections about the X — J plane. The lowest frequency bandgap 
for the even modes of the simulated cross structure extends 
from 0.91 GHz to 3.6 GHz. Within this bandgap, there are 
regions of full phononic bandgap (shaded blue) where no me- 
chanical modes of any symmetry exist, and regions of par- 
tial symmetry-dependent bandgap (shaded red) where out- 
of-plane flexural modes with odd symmetry about the x — j 
plane are allowed. As the bridge width is increased, the 
lower frequency effective-medium bands become stiffer, caus- 
ing an increase in their frequency, while the higher frequency 
tight-binding band frequencies remain essentially constant. A 
gap-map show in in Fig. [T]l, showing how the bandgaps in 
the structure change as a function of phononic crystal bridge 
width, illustrates this general feature. 

As shown in Figs. [2^-c, the cross crystal is used as a 
phononic cage (cavity) for an embedded optical nanocav- 
ity 1.3 3 J (highlighted in a green false color) consisting of 



a quasi-2D photonic crystal waveguide with a centralized 
"defect" region for localizing photons. This embedding of 
an optical cavity within an acoustic cavity enables, through 
the strong radiation-pressure-coupling of optical and acoustic 
waves, the probing of the properties of the bandgap-localized 
phonons via a light field sent through the optical nanocavity. 
The theoretical electric field mode profile and the measured 
high- 2 nature of the optical resonance of the photonic crys- 
tal cavity are shown in Figs.[2]l and [2^, respectively. Such a 
phonon-photon heterostructure design allows for completely 
independent tuning of the mechanical and optical properties of 
our system, and in what follows, we use this feature to probe 
arrays of structures with different geometric parameters. In 
particular, by varying the bridge width b of the outer phononic 
bandgap crystal, the lower bandgap edge can be swept in fre- 
quency and the resulting change in the lifetime, density of 
states, and localization of the trapped acoustic waves inter- 
acting with the central optical cavity can be monitored. Two 
different phonon cavity designs. Si and 5*2, were fabricated in 
this study (see Methods). We focus here on the lower acous- 
tic frequency Si structure, for which an array of devices with 
bridge width varying from b = 53 nm to 173 nm (in 6 nm 
increments) was created. Similar results for the ^2 structure 
are outlined in Appendix |B] 

Experimentally we observe the thermally excited acoustic 
modes of the photonic-phononic crystal through the induced 
phase-modulation of the optical cavity field |34|. The mixing 
of the phase-modulated light from the optical cavity with the 
transmitted light produces RF/microwave tones upon optical 
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FIG. 3: a, Optically transduced RF power spectral density of the 
thermal Brownian motion of Si structures with b = 57 nm (top), 
= 106 nm (middle), and Z? = 160 nm (bottom). The inferred re- 
gion below the phononic bandgap is shaded grey (see main text), 
b. Temperature dependence of the mechanical quality factor for the 
1.4 GHz acoustic mode of the Si structure with b = 51 nm. 



detection with a high speed photodetector (see Methods and 
Appendix |A]). The measured RF- spectra from three different 
structures with small (b = 57 nm), medium (b = 106 nm), 
and large (b = 160 nm) bridge widths are shown in Fig. [3^. 
Each narrow tone in the RF- spectra corresponds to a dif- 
ferent acoustic resonance interacting with the central optical 
nanocavity. Through careful calibration of the optical power 
and electronic detection, one can extract both the mechanical 
2-factor (from the linewidth) and the level of optomechani- 
cal coupling (from the magnitude of the transduced thermal 
motion) of each acoustic resonance. Here we parametrize the 
strength of the optomechanical coupling by the rate g, which 
corresponds physically to the shift in the optical cavity reso- 
nance frequency due to the zero-point motion of the acoustic 
resonance 1291 . A measurement of the temperature depen- 
dence of the acoustic mode spectrum is also performed, and 
is shown in Fig. [Sj) for one of the acoustic resonances over a 
temperature range from 300 K to 10 K. 

By measuring the entire set of Si devices in this way, a map 
may be produced of the localized acoustic modes' properties 
versus bridge width. In Figs. and |4}) we plot the numer- 
ically simulated and experimentally measured mode map for 
the Si structure. Each marker in these plots corresponds to a 
different acoustic resonance, with the position of the marker 
indicating the mode frequency and the size of the marker in- 
dicating the mechanical 2-factor of the mode (for the numeri- 
cal simulations all mechanical Q-f actors above 10^ are shown 
with the same marker size). Numerical simulations of the op- 



tical, mechanical, and optomechanical properties of the struc- 
ture are performed using the COMSOL [35 1 finite-element- 
method (FEM) software package, with an absorbing boundary 
condition applied at the exterior of the phononic cage | 27 , 3^. 
The various bandgap regions are indicated in Fig. with 
the same color coding as in Fig. [T]:. Due to the weak radi- 
ation pressure coupling to the optical nanocavity of the flex- 
ure acoustic modes of odd symmetry about the x-y plane of 
the slab (red mode bands in Fig. [T]:), we only show in the 
simulated mode plot of Fig. [3^ the even symmetry, in-plane 
acoustic resonances. 

The striking similarity of the simulated and measured mode 
plots is evidence that the optical nanocavity is able to sen- 
sitively probe the in-plane localized acoustic modes of the 
phononic bandgap structure (the acoustic band with light blue 
marker in Fig. is the one localized in-plane mode which 
does not show up in the measured plot of Fig. [4|); numerical 
simulations show this mode to be a surface mode at the inner 
edge of the cross crystal, which does not couple to the central 
optical cavity). Within the bandgap, modes are tightly local- 
ized (see Fig. [4]:) and do not radiate acoustic energy, whereas 
below the bandgap the acoustic modes spread into the exterior 
cross crystal (see Fig. [4]l), leaking energy into the surround- 
ing substrate region. The boundary where the mechanical Q- 
factor drops off is clearly identifiable in the experimentally 
measured mode plot of Fig.|4]3 (the spectral region below the 
apparent full phononic bandgap is shaded grey as a guide to 
the eye), and matches up well with the theoretical lower fre- 
quency band-edge of the full phononic bandgap of the cross 
crystal. 

Two other distinguishing features between modes inside 
and outside a bandgap are the spectral mode density and the 
strength of the optomechanical coupling. Below the phononic 
bandgap, acoustic modes can fill the entire volume of the cross 
crystal (out to the boundary of the undercut structure where it 
is finally clamped), resulting in an increase of the mode den- 
sity (proportional to volume) and a decrease in the optome- 
chanical coupling (proportional to the inverse- square-root of 
mode volume |29 1). In Fig. we have plotted the theoret- 
ically computed and experimentally measured values of the 
optomechancial coupling (g) for an acoustic resonance lying 
near the middle of the full phononic bandgap (this resonance 
is highlighted in orange in the theoretical plot of Fig. and 
red in the measured plot of Fig. The measured trend of 
optomechanical coupling nicely matches that of the theoret- 
ical one, and highlights the sharp drop off in optomechani- 
cal coupling as the mode crosses the bandgap. Similarly, in 
Fig. [4]f we plot the optomechanical coupling for the acoustic 
resonances of a single device with bridge width Z? ~ 100 nm 
(corresonding to a vertical slice in Figs.|4^ and |4j) as indicated 
by a dashed vertical line). Again we see good correspondence 
between theory and experiment, with the drop off in g and the 
large increase in spectral mode density clearly evident in both 
plots below the bandgap (note that the frequency position of 
the bandgap edge is not the same in the theoretical and exper- 
imental plots of Figs.|4^ and ^ due to the slight differences 
in bridge width). 

Having localized acoustic modes via a full 3D phononic 
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FIG. 4: a, Plot of the 3D-FEM simulated in-plane localized acoustic modes of the Si structure as a function of bridge width b. Each marker 
corresponds to a single acoustic mode, with the marker size proportional to the logarithm of the calculated acoustic radiation Q-factor. The light 
blue shaded markers correspond to acoustic bands which are optically dark. The shading corresponds to the same color coding of the phononic 
bandgaps as that used in Fig. [T}i. b, Measured mode plot of the optically-transduced localized acoustic modes for an array of Si structures 
with varying bridge width. The marker size of each resonance is related to the logarithm of the measured mechanical Q-factor. The inferred 
spectral region below the phononic bandgap is shaded grey, c and d, FEM simulations of the displacement field amplitude (|Q(r)|) for the 
acoustic mode in the orange colored band around 1.35 GHz in a. In c the mode is within the phononic band gap resulting in a radiation-limited 
Qm^^ ~ 10^. In d the mode is on the edge of the bandgap and has a reduced 2^^^^^ < 10^. e, Simulated (□) and measured (o) optomechanical 
coupling rate g for the orange (red) highlighted acoustic band in a (b). f, Simulated (□) and measured (o) optomechanical coupling rate g for 
the series of acoustic modes of the Si structure with b ^ 100 nm (vertical dashed curves in a and b) 



bandgap, and at least in principle removed radiation losses, 
it is interesting to consider the limits to mechanical damp- 
ing in these structures. As shown in Fig. [Sj) for the in- 
plane acoustic resonance at 1.40 GHz lying well within the 
phononic bandgap of the cross crystal, the mechanical Q- 
f actor increases from a value just above 3000 at room tem- 
perature to a value of 1.3 x 10"^ at a temperature of 10 K. 
This temperature dependence of mechanical Q is compared in 
Appendix [C] to Akheiser |37| damping, Landau-Rumer 1381 
damping, a numerical model of thermoelastic damping in the 
structure ||39l[39l|40l, and measurements of acoustic wave at- 
tenuation in bulk Si ll4Tl . The measured temperature depen- 
dence and the overall magnitude of the measured mechanical 
2-factor are seen to be much smaller than any of these com- 
parisons, suggesting that surface effects and/or fabrication- 
induced damage may be playing an important role in the me- 
chanical damping of the nanosctructed devices studied here. 

Beyond the confinement and localization of acoustic modes 
in three dimensions, the connected geometry of the 2D-0MC 
structures presented in this work offers a platform for more 
complex phonon-photon circuitry. As has been described in 
recent theoretical analyses l29l|3T]|32J, such circuitry could be 



used to create optomechanical systems with greatly enhanced 
optomechanical coupling, and to realize devices such as trav- 
eling wave phonon-photon translators and slow light waveg- 
uides f42l capable of advanced classical and quantum optical 
signal processing. The functionality of these devices are based 
upon the slow propagation velocity and long relative lifetime 
of phonons in comparison to photons, which allows for the 
storage, buffering, and narrowband filtering of optical signals. 
In addition, the coupling of optomechanical circuits to a wide 
variety other physical systems, such as superconducting elec- 
tronic circuits B3l and atomic vapors ll44ll . may also enable 
the interfacing and networking of different quantum systems. 
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Methods 

Fabrication: Phononic-photonic cavities are fabricated using 
a Silicon-On-Insulator wafer from SOITEC (p = 4-20 ^1-cm, 
device layer thickness t = 220 nm, buried-oxide layer thick- 
ness 2 /im). The cavity geometry is defined by electron 
beam lithography followed by inductively-coupled-plasma re- 
active ion etching (ICP-RIE) to transfer the pattern through 
the 220 nm silicon device layer. The cavities are then under- 
cut using HF:H20 solution to remove the buried oxide layer, 
and cleaned using a piranha/HF cycle f45\. 
Optical measurement technique: Optically we character- 
ize the samples using a dimpled fiber taper setup connected 
to a swept- wavelength external-cavity laser |46|. A typical 
fiber taper transmission spectrum is shown in Fig.[2je), with a 
measured intrinsic optical quality factor of Qi = 1.5 x 10^. 
By controlling the taper position when touching the sam- 
ple we are able to control the coupling between our fiber 



taper-probe and the optical cavity. Usually, after touching, 
the external coupling rate was on the order of tens of MHz 
(Ke/27r ^ 70 MHz) which corresponds to a transmission dip 
of ^70%. 

RF measurement technique: The transmitted cavity light is 
sent into an erbium doped fiber amplifier (EDFA) and then 
onto a high-speed photodetector. The photodetected signal is 
sent to an oscilloscope (2 GHz bandwidth) where the elec- 
tronic power spectral density (PSD) is computed. Since our 
devices are in the sideband resolved regime, i.e., the total op- 
tical loss rate, /c, is smaller than the mechanical frequency, 
Q,M, the largest transduced signal is achieved when the laser 
frequency is detuned from the optical cavity resonance by ap- 
proximately the mechanical frequency |11 1. Therefore, for all 
the measurement, the probe laser was locked ^ 1 GHz away 
from the cavity resonance for Si structures and ^2.5 GHz for 
5*2 structures, both on the blue (higher frequency) side of the 
cavity. 
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Appendix A: Experimental Setup 

The experimental setup used to measure the phononic- 
photonic crystal cavity properties is shown in Fig. [5ja). A 
fiber-coupled tunable infrared laser, (New Focus Velocity, 
model TLB -6328) spanning approximately 60 nm, centered 
around 1540 nm, has its intensity and polarization controlled 
respectively by a variable optical attenuator (VOA) and a fiber 
polarization controller (FPC). The laser light is coupled to a 
tapered, dimpled optical fiber (Taper) which has its position 
controlled with nanometer- scale precision. The transmission 
from the fiber is passed through another VOA before being 
detected. 

To measure the optical properties, a photodetector (PD, 
New Focus Nanosecond Photodetector, model 1623) is used. 
The detected optical transmission signal is recorded while 
sweeping the laser frequency. By controlling the distance be- 
tween the fiber taper and the sample, the external coupling 
rate (/q) is changed. Fig. [5jc) shows the change in the cou- 
pling rate for two different positions of the fiber taper. In the 
limit where the external coupling rate is zero we can measure 
the intrinsic coupling rate (Ki). The total optical loss is then 

K= Ki. 

To measure the mechanical properties, the transmitted sig- 
nal is sent through an erbium doped fiber amplifier (EDFA) 
and sent to a high-speed photoreceiver (PR, New Focus 
model, 1554-B) with a maximum transimpedance gain of 
1,000 V/A and a bandwidth (3 dB rolloff point) of 12 GHz. 
The RF voltage from the photoreceiver is connected to the 
50 input impedance of the oscilloscope. The oscilloscope 
can perform a Fast Fourier Transform (FFT) to yield the RF 
power spectral density (RF PSD). The RF PSD is calibrated 
using a frequency generator that outputs a variable frequency 



sinusoid with known power. 

As stated in the main text, our devices are in the sideband 
resolved limit, i.e. the total optical loss rate is smaller than 
the mechanical frequency, K < Q^u- Therefore the largest 
transduced signal is achieved when the laser frequency de- 
tuned from the optical cavity resonance by the mechanical 
frequency fTTl. The probe laser is locked to approximately 
1 GHz (2.5 GHz) on the blue side of the cavity resonance for 
the Si (S2) structures. By measuring the transmission contrast 
during the acquisition of the RF PSD and comparing with the 
transmission curve of each device (as shown in Fig. [JJc)) we 
determine the laser detuning and the dropped power into the 
cavity. To lock the probe laser frequency a given frequency 
away from the resonance, a 90/10 beam splitter (BS) is added 
to the optical path, and the signal from the 10% arm is feed 
to a PD connected to a locking circuit which compares the 
voltage level from the transmission signal to a predetermined 
value to generate an error signal. The error signal is then fed 
into the laser to stabilize the laser frequency. 



Appendix B: Phononic Band Gap Tuning via Lattice Size 

In the main text we demonstrate a phononic bandgap 
around 1.3 GHz based upon the Si design. In this section we 
present measurements performed on a second design (5*2) with 
phononic bandgap in a different frequency range. Following 
the nomenclature in Fig. 1 of the main text, the nominal di- 
mensions for the 5*2 devices are: a = 925 nm, h = 850 nm, and 
w = 210 nm, which allows for an acoustic bandgap around 1.6 
to 2.5 GHz. The optical nanocavity in the 5*2 devices are the 
same as that used in Si devices. 

To show the presence of a bandgap, we fabricate and mea- 
sure a set of devices with bridge widths varying from 70 to 
170 nm. The RF PSD of the optically transmitted signal is 
shown in Fig. [6^ for three different ^"2 structures with different 
bridge widths. By comparing the spectra for different bridge 
width the presence of a phononic bandgap becomes evident. 
Figs. [6]3 and [6]: shows the simulated and measured acoustic 
mode plots for the in-plane acoustic modes of the 5*2 set of 
devices as a function of the bridge size b. As in the main 
text, the marker size is proportional to the mechanical qual- 
ity factor. We can easily identify the lower bound of the full 
bandgap region (gray shaded area) by the abrupt change in the 
mechanical quality factor and the mechanical density of states. 
Similar to the Si devices, as the bridge width is increased, the 
lower frequency bandgap boundary tunes to higher frequency. 
The mechanical mode localization is also affected by the pres- 
ence of the bandgap as seen in Fig.[6]i and Fig. [6^. There the 
displacement field amplitudes of a confined mode are plotted 
within (Fig.|6]l) and on the edge of the bandgap (Fig.|6^). 

Appendix C: Mechanical Damping 

For micro- and nano-mechanical structures the many 
loss mechanism include clamping losses, dissipation due 
to phonon-phonon interactions, and surface effects. In the 
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FIG. 5: Experimental Setup a, Experimental setup for measuring the PSD. b, Optical micrograph of the tapered fiber coupled to the device 
while performing experiments, c, Optical spectra for two different positions of the taper relative to the device are shown. 



present work, the engineered structures virtually eliminate the 
clamping loss for modes within the full mechanical bandgap. 
As a result other loss mechanisms can be studied. In the 
main text, we present a series of measurement of a con- 
fined mode, at Q^u/'^^ =1.4 GHz, versus the experimental 
ambient temperature which shows a temperature dependent 
mechanical quality factor, 2m. In Figure |7ja) we compare 
our results (blue circles) with other temperature-dependent 
losses. For the acoustic attenuation, a{Q,',T), the relevant 
bulk measurement results are taken from Ref. | 41 1 and we plot 
(black diamonds) the calculated Qu = 27iQ.m/{'^oc), where 
^Im/^tt =1.4 GHz is the mechanical frequency. The mea- 
sured values from Ref. |41 1 indicate that our devices are not 
limited by the bulk losses. 

In order to provide an upper bound for the quality fac- 
tor in our structures, we compare our measurements with the 
temperature-dependent acoustic attenuation provided by the 
Akheiser 1 37 1 (green dashed line), Landau-Rumer 138 ! (pur- 
ple dashed line), and Thermoelastic Damping (391 SOI 133 
(TED - red square points) models. The Landau-Rumer model 
provides a microscopic theory for sound absorption and is 
valid on the limit of £1m^ ^ 1, where T is the mean time be- 
tween collisions of thermal-phonon in the solid. Akheiser' s 
model treats the dissipation of heat generated by strain trough 
the Boltzmann transport equation [48 1 and is valid for Q^u^ <C 
1 . Finally TED considers the heat diffusion in homogeneous 
materials which can be calculated for any geometry. 

The equations for acoustic attenuation for Landau-Rumer 
(oclr) and Akheiser's (aAx) models used are |48|: 



and ap,Y.(Q.\T) = 



(Cl) 



where y is the average Griineisen coefficient extract from 
Ref. ll49l , Cv is the volumetric heat capacity with values from 
Ref. |W|, and p = 2330 kg/m^ and q = 9.15 x 10^ m/s are 
the density and speed of sound of Si respectively. Fig. [TJa) 
shows that for temperatures above T = 100 K the Akheiser 
and Landau-Rumer losses are dominant. Note that Q.m^ = 1 
for T = 200 K and only Landau-Rumer is valid below this 
point. In this region TED effects due to the sample geometry 
become important. 



Our approach for calculating TED follows that of Ref. fSTIl. 
where the TED-limited 2m,ted is extracted from 2D-FEM 
simulations 1 351 for the thermo-mechanical equations consid- 
ering a finite thickness. Fig.|7jb) shows the change in temper- 
ature, Ar(r) = r — r^, for the deformed structure, versus the 
phase of the mechanical oscillations, at the expected brownian 
motion amplitude for = 300 K. The temperature difference 
was calculated based on the maximum thermal displacement 

amplitude Xmax = y^^H^/m^ff^l^. These plots show that dur- 
ing a mechanical cycle, even when there is no deformation, 
i.e. ^^M^ = 7i/2 and 3;r/2, the temperature gradient is non- 
zero. This indicates that the temperature does not follow adia- 
batically the strain/stress profile, causing a time-delayed force 
to be imparted on the resonator, and leading to dissipation. 
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FIG. 6: S2 cavity design, a, RF-power spectral density for three devices with different bridge widths, h. The gray shaded area refers to 
regions outside of the full bandgap while the orange shaded area are for frequencies within the bandgap. b, results of full 3D-FEM simulations 
of phononic localized modes as a function of bridge width h. For each bridge- width a full 3D-FEM simulation with an absorbing perfect- 
matched-layer (PML) is performed. Each square corresponds to a single mode, with its size is proportional to the logarithm of the mechanical 
2-factor. The gray shaded squares represent modes which are optically dark (g's small enough to be lower than our detection noise floor), c. 
Measurement (o) results of the localized mechanical modes for the a set of fabricated ^2 samples. Each bridge size h was measured by careful 
analysis of scanning electron micrographs. The marker sizes are proportional to measured mechanical Q-f actor, d and e, FEM simulations 
of the displacement field amplitude |Q(r)| for the mechanical mode at 2.1GHz shown in b. In d, the mode is within the phononic bandgap 
resulting in a radiation limited 2^^^^^ ~ 10^. When on the edge of the bandgap, the mode is less localized as shown in e, and has a reduced 
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FIG. 7: Temperature-dependent mechanical losses, a, Comparison between the different sources of mechanical loss with the measured 
Qm values versus temperature. The circles represent the measured values from the Si samples; diamonds are computed 2m from acoustic 
attenuation data from Ref. |41 1; squares represent the simulated values for TED as explained on the text. The purple and green lines are the 
calculated Qm due to Akheiser and Landau-Rumer phonon-phonon dissipation mechanism respectively. The insets show the measured RF 
PSD at 10 K and 300 K for extracting Qm- b Thermo-mechanical 2D-FEM simulations for the mechanical mode at 1.4 GHz shown in Fig. 3 
of the main text. The thermal profile is plotted at various times during the mechanical cycle. 



